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A geometrical model for tri- nuclear molecules is presented. An analytical solution is 
obtained provided the nuclei, which are taken to be prolately deformed, are connected 
in line to each other. Furthermore, the tri-nuclear molecule is composed of two heavy 
and one light cluster, the later sandwiched between the two heavy clusters. A basis is 
constructed in which Hamiltonians of more general configurations can be diagonalized. 
In the calculation of the interaction between the clusters higher multipole deformations 
are taken into account, including the hexadecupole one. A repulsive nuclear core is 
introduced in the potential in order to insure a quasi-stable configuration of the system. 
The model is applied to three nuclear molecules, namely 96 Sr + 10 Be + 146 Ba, 108 Mo + 

10 Be + 134 Xe and 112 Ru + 10 Be + 130 Sn _ 



X 



1. Introduction 

In a recent experiment Q the cold fission decay of 252 Cf into three clusters was inves- 
tigated. In coincidence measurements the three participants were identified as being 
96 Sr, 10 Be and 146 Ba. In reanalyzing the data, two further systems where discovered 
namely 90 Y + 142 Cs and 108 Mo + 134 Te, all with 10 Be as the third particle. Surprisingly 
the data suggest that the transition from the first excited 2 + state to the ground state 
in 10 Be was not Doppler broadened as one would expect if the system immediately sepa- 
rates into three clusters and the Be nucleus deexcites in flight. In addition one observes 
that the transition energy of 3368 keV in the 10 Be nucleus interacting with 96 Sr and 
146 Ba is probably by about 6 keV smaller than for the free 10 Be nucleus. The transition 
energy decreases further for the other two systems, being largest when both heavy clus- 
ters are spherical. The heuristic explanation is that the average shell model frequency 
in presence of the two heavy clusters is modified. The influence of both clusters leads 
to a softening of the 10 Be potential and thus to a somewhat smaller transition energy. 
The largest overlap, i.e. the strongest change in the average shell model frequency, of 
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one heavy cluster with Be is obtained for a spherical deformation. Though, one has 
to wait for further experimental confirmation and better statistics, the indication of 
the possible existence of three-cluster molecules calls for a theoretical description and 
prediction of their structure. 

The interpretation of the above observation is the probable existence of a nuclear 
molecule ||] with a half life larger than 10 -13 sec |IJ. Such large lifetimes would open 
up the possibility of a spectroscopy of giant nuclear molecules. In Q we proposed a 
phenomenological, geometrical model for the system of three clusters, two heavy ones 
and a light in the middle. The model was restricted to butterfly and belly-dancer modes 
only and the 10 Be cluster was assumed to be spherical. However, in general the light 
cluster can be deformed and its effect must be studied. Also the inclusion of (5 and 
7 vibrations must be considered. Because the light cluster in |3j was considered to be 
spherical, the stiffness of the butterfly motion is mainly determined by the monopole 
part of the Coulomb repulsion between the heavy fragments []|]. In the case when the 
light fragment is deformed, this is no longer true and one has to determine explicitely 
the change of the nuclear and Coulomb interaction between the light and heavy clusters 
as a function of the inclination angle. The effective potential is determined via a double 
folding calculation and is described in Ref. ||]. In Ref. ||] the nuclear potential was 
taken into account, including also multipoles higher than the monopole and quadrupole 
ones. 

In this contribution we present the details of the model proposed in Ref. Q which 
itself is an extension of Ref. presented for two clusters. The molecule exhibits 
butterfly and belly-dancer modes, ft- and 7-vibrations of the two clusters. This picture 
can be extended straightforwardly to three clusters using the formulas as given in Ref. 
Q. It is especially easy for the case when the two big clusters (for example, the 96 Sr and 
146 Ba) are connected via a smaller spherical nucleus ( 10 Be). The situation is illustrated 
in Fig. 1, where the main dynamical variables, with a spherical cluster in the middle, 
are indicated. In this paper we will extend the model to include also (3 and 7 vibrations, 
with the addition of a deformed light cluster. 

In Fig. 2 possible vibrational modes are indicated, when the light cluster is deformed. 
The first case corresponds to the above mentioned butterfly mode, whereas the second 
one is the antibutterfly mode and is not investigated in this contribution. We will see 
that for the butterfly motion the dominant contribution comes from the movement of 
the center of masses and, contrary to the two cluster case, deformations play a minor 
role, except in fixing the length of the axis in Fig. 1 and 2. Also the methods presented 
in this paper can be readily extended to this antibutterfly mode. 

In order to keep the problem tractable, the main assumptions are that the light 
cluster is sandwiched in-between the heavy nuclei and the inclination angles of the 
clusters are small with respect to the molecular (fission) z-axis, connecting both heavy 
clusters. In such a linear chain configuration the total potential of the lighter cluster 
has an absolute minimum on the axis joining the three fragments. As can be seen in 
Fig. 3, for a given distance d between the tips of the two heavier fragments, there is a 
point on the fission axis z, where the forces exerted by the heavy fragment 1 on the light 
cluster are canceled by the forces exerted by the heavy fragment 2. This is the so-called 
electro-nuclear saddle point The details of the potential calculation are explained 
in section 3. For the moment we just note that the result of Fig. 3 was obtained using 
a strong nuclear repulsive force, between the fragments in order to avoid their mutual 
overlap. 

The paper is organized as follows: In section 2 the general formulation of a three 
cluster molecule is presented. For practical purposes we then restrict to the case where 
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the three clusters are in touch and all three can be deformed. This section will be 
divided into one where no (3 and 7 vibration is included and a second one where these 
degrees of freedom are considered, too. In section 3 the derivation of the nuclear and 
Coulomb potentials is given. In section 4 the model is applied to the three systems: 96 Sr 
+ 10 Be + 146 Ba, 108 Mo + 10 Be + 134 Te and 112 Ru + 10 Be + 130 Sn. All these splittings 
are a product of the ternary cold fission of 252 Cf. The structure of the participants is 
discussed. Finally, in section 5 conclusions are drawn. 



2. The general formalism 

The motion of the three clusters can be divided into the rotations of the individual 
clusters plus the motion of their center of masses with respect to each other. The part 
of the Hamiltonian which describes the individual rotations can be read off from Eqs. 
(49) and (50) of Ref. 0. The part which changes originates from the motion of the 
center of masses. Therefore, the discussion concentrates first on the motion of the three 
center of masses and is independent of considering either of the two modes described in 
Fig. 2. 

In order to separate the center-of-mass motion, the following coordinates are intro- 
duced 



r 2 — ri 
m\r\ + m 2 r 2 

mi + 1^2 



rs (1) 



where = A^m, is the number of nucleons of cluster no. k and m the nucleon 
mass. The first coordinate (r) describes the relative distance of the two heavy fragments 
while the second one (£) the distance of the third, lightest cluster to the center of mass 
of the first two. The kinetic energy, excluding the motion of the total center of mass, 
acquires the form 

T = + ^(12)3^ 2 + l^UlJlUl + \ t U 2 J2U 2 + ^U 3 J 3 UJ 3 (2) 

where H12 = 2 (ml+^ 2 ) and ^(12)3 = 2 (Z ( +mt+m 3 ) ■ The first term in @ describes the 
kinetic energy of the two heavy clusters with respect to each other and the second term 
the kinetic energy of the third cluster with respect to a mass (mi + m 2 ) at the center 
of mass of the first two clusters. The mass factors describe the reduced mass for each 
case. The term proportional to f 2 has the same form as for the two cluster case and 
thus is already included in the considerations of Ref. . The last three terms in eq. @ 
are describing the rotational motion of the three clusters with angular velocities ^1,2,3, 
referred to the laboratory frame. The inertia tensors Si are defined in the intrinsic frame 
such that in the absence of (5 and 7 vibrations the only non-vanishing components are 
the first two diagonal terms, {JJi)\\ = (J~i) 22 = Ji, the quantum rotation around the 
symmetry axis of any of the two heavier fragments being discarded. When /3 and 7 
vibrations are included there will be a contribution to (^4)33 given by a 7 dependence 

The second term in Eq. @ needs more attention. In Fig. 4 the three center of masses 
are plotted and the relevant coordinates are indicated. The projection of the vector £ 
onto the relative distance vector, denoted by £ 2 , and its perpendicular component along 
the x-axis (£ x ) are given by 
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(,x = n3sine (3) 

where r%z is the distance between the cluster no. 1 and 3 and e is the angle between the 
axis connecting cluster 1 and 3 to the vector r. Note, that this angle is not necessarily 
the same as the inclination angle of the intrinsic z-axis of cluster no. 1 to the vector r. 
Because the molecular plane is defined by the three center of masses of the clusters the 
spherical components £±i are the same up to a sign (see Appendix A and note the change 
in sign in the definition of the spherical components, due to convenience, compared to the 
usual definition H ) . The £ contribution of the kinetic energy is obtained by rotating into 
the molecular system and then substituting the expressions for the spherical components 
of the vector In order to keep the procedure manageable, we assumed that the angle 
e is small, i.e. the light, third cluster is not far away from the axis connecting the 
two heavy clusters. Furthermore, we assume that the first cluster, supposed to be the 
lightest one of the two heavy clusters, touches the third one and the third cluster touches 
the second one. The situation is illustrated in Fig. 5 with a certain exaggeration as 
concerned to the distance of the third cluster to the axis (r) connecting the heavy 
clusters. These assumptions exclude the anti butterfly mode for which the relation of 
the angles change. The procedure for that mode, however, is completely analogous. 

In Fig. 5 the relation of the angle (p 3 to p\ and p 2 can be read off with the supposition 
that the three clusters are connected. Assuming small angles, the result is 

V? 3 ps sin p 3 = — - (R 2 sin ip 2 - Ri sin <p x ) rj — — (R 2 ip 2 - Ri <Pi ) (4) 
ZK 3 ZK 3 

In the above formula we suppose that the heavy cluster 2 has a larger shift on the in- 
direction than the lighter one. In principle the two angles ipi and (p 2 could be treated 
now as independent with <p 3 constrained by Eq. (|4]). However, the problem would get 
too complicated, implying coupling terms between the ipi and tp 2 motion. Supposing 
that the cluster in the middle is small, the relation between <p\ and ip 2 should not differ 
much from the case when the small cluster is spherical We use, therefore, the same 
relation ip 2 m^rnf] <Pl as f° r a spherical cluster and substitute it into Eq. (||) resulting 
in (ipi e) 

One possibility to relax this constraint is to expand the general motion around the 
point where condition (||) holds and diagonalizes the Hamiltonian in this basis. This 
standpoint will be adopted in the present paper. 

For the case of connected clusters the distance between center of masses of the 
third cluster and nucleus no.l, R13, can be approximated by (R\ + R3) in the linear, 
unperturbed configuration. In the perturbed case, when the molecule ceases to be linear, 
rx3 is modified by the amount 

1 R1R3 {Ri + R 2 + 2R 3 \ 2 2 

5ri3 * ~8r7Tr3~ { WTr, ) £ (6) 

provided e is small, Similarly we get the variation of r 23 



1 R 2 R 3 fRi + R 2 + 2R 3 \ 
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5r ™ W SR^TW, { R 2 + R 3 ) £ (7) 
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At equilibrium (e=0), the distance r, between the center of masses of the two heavy 
clusters is given by Rq = R\ + R% + 2i?3. When £ / 0, and taking into account terms 
up to second order in e, r changes to 

r^(R 1 + 2R 3 + R 2 ) (l - I|±±i| e 2 ) + 6r 13 + 5r 23 . (8) 

Therefore, an increase in e, changes r, ri3 and r23 with a correction in the moment 
of inertia of the order of e 2 and higher. At their turn, the components £ z and £ x acquire 
the form 



!l> (R 1 + 2i? 3 + R 2 ) - (Ri + R 3 ) 



(A 1 + A 2 ) 



^ « (fli + fla)e + *Cx ■ (9) 

where <5£ 2 and (5^ give the contributions due to changes in the relative distance of the 
clusters 1 to 3 and 2 to 3, i.e. they describe a mode of the stretching vibrations. In the 
unperturbed configuration 

. A 2 (R 2 + R 3 )-A 1 (R 1 +R 3 ) 

e ° = mTm (10) 

Again, for simplicity we assume that the relative vibrations are along the molecular z 
axis, i.e. 5£ x = 0. If this assumption is not made, there will be contributions of the 
type 5^ x 5^ z in the kinetic energy. One has then to construct the basis first, where this 
coupling term is absent, and afterwards to diagonalize the complete Hamiltonian in 
this basis. Note that the model is quite crude and other coupling terms, like rotation- 
vibration interactions, have been neglected. Consequently, the spectrum obtained in the 
model will be a first approximation and the assumptions made above are justified only 
in this context. 

The radial mode £ along the z axis describes the motion of the two heavy clusters 
in a common direction and of the light cluster in the opposite direction (see Ref. Q). 
The other radial mode comes from changes in r and describes a vibration of the two 
heavy clusters with respect to each other while the small cluster does not move. These 
types of vibrations will be included in the complete Hamiltonian. Inserting (^) into the 
kinetic energy for the £ motion, using 5^ x = 0, we obtain(see Appendix A) 



Tt = \ii {l2)3 {(Ri + R 3 fe 2 + e^? + ^ 2 2 ) + {Ri + Rz) 2 e 2 J 3 



I /.,■■] .-)-,■) < 1 n 1 -- 1— n •> 1 -- 1— r,,\ ;>;■■■. -i- i -i- \ n i -\- n •..> i .■ 

2^(R 1 + R a )(iu' 2 -eu' 1 ui i )} , (11) 



with oj'u being the angular velocity around the fc'th molecular axis (1 = x, 2 = y and 
3 = z) and £ is a shorthand notation for 5^ z . From here on, when we refer to the total 
contribution, including the rotational one, we continue to denote it by £ and the pure 
motion along the z axis is denoted by £. 

The kinetic energy of the r and £ stretching vibrations along the molecular z axis is 
described by 

T r = \llX2T 2 + ^(12)3^ , (12) 
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where the last term comes from Eq. (p"l|). The coordinates r and £ are related to the 
relative distances 7*13 = (213 — z\) and r23 = (z 2 — Z3) via 



(Z2 ~ Z 3 ) = 
(>3 - Zl) = 

The corresponding potential is given by 

V r = -7p(( r 3 - n) - ri 3i0 ) 2 + -|^(( r 2 - ^3) - r 2 3,o) 2 (14) 

where ri3,o and r23,o are the equilibrium positions 7*13,0 = (R1+R3) and 7-23,0 = (-R2+-R3) 
respectively. Note that in our picture cluster no. 1 is to the left and no. 2 to the right 
and that the coordinates refer to the distance along the intermolecular axis. 

Up to know we considered the r and the £ motion only. The contributions coming 
from the deformation of the individual clusters can be read off from Ref. and from 
there the general kinetic energy can be constructed. 

In the next subsection we will consider two different cases mentioned in the intro- 
duction. In the first one the contributions of (3 and 7 vibrations are excluded and in 
the second one they are included. In both cases the static deformation in the ground 
state is assumed to be prolate. This restriction is governed by the necessity to keep the 
problem solvable otherwise the complicated form of the Hamiltonian would prevent an 
analytical solution. In case a triaxial nucleus is present the procedure outlined is strictly 
speaking not valid, but an approximation of the triaxial nucleus by an axial symmetric 
one would do the job. 

2.1 Without (3 and 7 vibrations 

With the assumption made in the preceeding section the total kinetic energy is given by 



mi 



m i + m 2 
m 2 

mi + m 2 



r + £ 
r — £ 



(13) 



T = -OllK 2 + J 2 2 ) + ^633^3 2 " 9 W<"3 + \®eS 2 + &2siu' 2 
1 1 

with 



8n Ri Jl + J 2 + h+ M(12)3^0 + Ml2-Ro 

Q 2e = 9i 3 (16) 

For nearly symmetric heavy clusters, not too deformed and a small value of m ™^ m2 the 
coupling is small compared to the diagonal terms of the moment of inertia and thus 



6 



can be neglected, as for the other contributions (see Ref. (Tj). Also the deformation 
dependent part in 6n can be neglected with respect to the last term, as was done also 
in Ref. 0. However, for very asymmetric heavy clusters we cannot neglect any more 
the contribution of 02 £ , except for small m ^+ m ■ in this case, one has to diagonalize 
this term in the basis with ®2e = 0, which is in the same spirit as for the coupling terms 
in the radial vibrations. The basis will be discussed further below. 

For the potential we assume a quadratic behaviour in e, r and £, i.e. 



V = ^ 2 + §- 2 + ^-£o) 2 , (17) 

where f = (r — tq) with ro being the equilibrium position of the nuclear molecule. The 
parameters C r and are related to C13 and C23 via 

c _ (m 2 2 C 13 + mf C 23 ) 
(mi + m 2 ) 2 

= C13 + C23 . (18) 

Other crossing terms of the type r£, r and £ also appear, which vanish for a symmetric 
dinuclear sub-system, formed from clusters 1 and 2. We assume that in general the 
microscopic interaction is such that also for non-symmetric clusters the coupling terms 
can be canceled or made small, which is obviously a simplification. 



In what follows, we quantize the Hamiltonian with the kinetic energy given in (15) 



taking into account the contribution of the coupling of the e and u>2 motion, whose origin 



is the Coriolis force. The quantization procedure is explained in Ref. [10] and was also 
used in Ref. @. 

Then a metric can be identified via 

2Tdt 2 =Y,9^dq^dq u (19) 

where T is the total kinetic energy and dq^ is a short hand notation for the variables 
appearing in our model. The quantized kinetic energy is given by 

= -Y^W^ ^ (20> 

with g = det(g fJiU ). The volume element is 



dr = yfgdqi ... dq n = ^fgdr' . (21) 

It is more convenient to use as volume element dr' instead. We have to redefine the 
Hamiltonian H and the wave function of the original Schrodinger equation via 

H -» H' = g^Hg~^ 

tfj -> <f> = gh (22) 

so that 



J ^*Hi)dT = J cp*H'cpdT' 
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holds. The potential commutes with g because it does not contain differential operators. 
The kinetic energy simplifies to 



h 2 t-^ d d 

>^dq~ u + Vadd 



V, 



add 




1 dg dg 1 d{g 1 



^ g 2 dq u dq v g dq^ dq v 



g dq u dq 



The Hamiltonian composed by the kinetic energy (|T^) and the potential is 
quantized in this manner, resulting for the kinetic energy in 




+ TS^TT T^-S4 7Z (24) 



In the following discussion we will skip the term in the parenthesis {...} which has to be 

02 

treated as a perturbative interaction and the corrections given by e ^ , which is small 
compared to one for the molecular systems studied in this paper. In very assymmetric 
systems, however, both terms have to be included and diagonalized in the basis which 
will be constructed in what follows. The terms which do not contain a derivative come 



from the additional potential of Eq. (23) 



Neglecting the term with the parenthesis {...} in Eq. ( |24| ) and corrections of the 

Q2 @2 

order of q^- and -g^-, the corresponding static Schrodinger equation can be solved with 
the ansatz 

4> = DMK^)XK,nM9n r (r) (25) 

were 

d mk($) is the Wi 

gner D-matrix g nr ( r ) is the one dimensional harmonic oscil- 
lator for the relative motion and XK,n E is the solution of the differential equation 

h 2 d , J_( K 2 _ i\h 2 , c^ji 
20 e dp " r £ ^ 4 ' 2~e^ 2 fc 

+ ( K ^ i [I{I + l)-K 2 ] + hio r {n r + \) + huj^n l : + } 1 )-E^\ <j> =0 , (26) 



where I is the total spin, K its projection to the molecular z axis, huj r = foy ^~ > 

= h\ ,, * and E the total energy of the state 6. 

y M3(i2) aj r 

This equation can be solved with the solution given by (see also Ref. W~ 



{\ lk+ 2T(l K + ^ + n £ )) 2 i, a 3 

X^(e) = ] - - J ^- +1 e-2^ 1 F 1 (-n e ,ij f + -;Ae 2 ) , (27) 

(n E !)sr(Zjr + |) 2 
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where A 2 = I 



K 



K 



and \Fi{...) is the confluent hypergeometric function. 



The total energy is given by 



E 



26 
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I(I + 1)-K 2 +hu r (n, 



+ l) + hu 6 QK\+2n e + ^) (28) 



and hco e = h\ 

When there are two identical clusters the wavefunction has still to be symmetrized 
with respect the interchange of cluster no. 1 and 2. For details on how the variables 
change see Ref. at least for the two heavy clusters. For the light cluster the variables 
transform in the same way as indicated by the second cluster except for the change of 
indices from 2 to 3. With respect to the variable e the transformation is e — > —e. 



2.2 Including j3 and 7 vibrations 

The extension to (3 and 7 vibrations is straightforward. Including also the rotation 
around the intrinsic z axis of axial symmetric nuclei, given by (k = 1, 2, 3), leads to 
a very complicated form. Additionally we have to assume that the clusters are prolately 
deformed. Otherwise a complex coupling between the rotations around the x, y and z 
will appear. The (3 and 7 vibrational variables of the fe'th cluster are defined by 



Pk — a o — Pok 

Vk = 4 , (29) 

where Oq and are the components of the quadrupole deformation variable of the fc'th 
cluster with respect to the principal axes. Using Eq. (50) of the second paper in Ref. 
Q for the moments of inertia plus the corrections discussed further up, the Hamiltonian 
acquires the form 




lGBir/l 16B 3 r% IQB^rjl 

n 2 d 2 h 2 d 2 h 2 d 2 



4Bi d 2 r,l 4B 2 d 2 r)\ AB 3 d 2 rg 

h 2 d 2 h 2 d 2 n 2 d 2 



lQB ir j 2 d<S> 2 IQB2V2 d ®l !6^3??3 <9$| 

h 2 d 2 h 2 d 2 

2/i(i2) dr 2 2/i (12)3 d£ 2 



(30) 
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where the variables fa and rjk (k = 1,2,3) describe the (5 and 7 degree of freedom as 
defined in Ref. @. 

The complete potential is given by 



v = ^ e 2 + ^f 2 + ^%-£ ) 2 

2 2 ^ 2 K ' 



+C m rj 1 + C. m r\ 2 + C V3 rj 3 

1 52 , C/32 52 1 C/33 52 



# + ^2 + ^/9g (31) 



and the determinant g is 



5 = S^Q^eiBtBiBirjiv^ie 2 . (32) 

where 9 22 ~ ©n was used. The factor 2 3 B 2 B%B 2 = (2Bi)(2B 2 )(2B 3 ) BiB 2 B 3 comes 
from the (3 and 7 part of the metric tensor g^ u . 

Neglecting, as in the case without (3 and 7 vibrations, corrections of the order of 

02 

q e q x1 the static Schrodinger equation can be solved with the ansatz 

9n ( (O9n 01 (fa)9n 02 (fa)9ng 3 (PtiXKi,^ {Vl)XK 2 ,n V2 (m)XK 2 ,n V2 M (33) 

where is the eigenvalue of the operator ijS— and K stands for I K — K\ — K 2 — K 3 I. 
Applying the transformations of the coordinate symmetries, as explained in Ref. M and 
Appendix B, the final form of the wave function, given in Appendix B, is obtained. 
The energy is 

h 2 

E = -—{I(I + l)-K 2 ]+Tiu Je {\K-K 1 -K 2 -Kz\+2n e + l) 

111 

+hu Pl (n Pl + -) + hujp 2 {n{3 2 + -) + hufo{np 3 + -) 

+hu Vl (- | K\ | +2n Vl + 1) + fiu^ m {^ I K 2 \ +2n m + 1) 
+hu V2 (~ | K 3 | +2^3 + 1) + ?iu r (n r + -) + fiu)z(ri( + -) . (34) 
The frequencies are 



hu r = hj-p^, fuj* = hJJ!s- . (35) 

Before we apply the outlined procedure to the cases 96 Sr + 10 Be + 146 Be, 112 Ru + 
10 Be + 130 Ru and 108 Mo + 10 Be + 134 Te, which are all non-symmetric systems, the 
derivation of the potential parameters is outlined. The ones related to the [3 and 7 
vibrations are obtained through the spectrum of the individual clusters. For details see 
Refs. 0, [L0|. 
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3. Derivation of relative potentials 



The aim of this section is to derive the expressions for the stiffness coefficients appearing 
in Eq. dl7|) . For that we need to calculate the interactions between the nuclei composing 
the trinuclear molecule. Since the interaction potential between the clusters should 
depend not only on their reciprocal distances but also on orientations we choose a 
double folding-model potential in which the nuclear densities p are directly folded with 
effective nucleon-nucleon interactions : 



V(R) 



dr^ 



dr 2 pi(ri)p 2 (r 2 )v(s) 



(36) 



For the nuclear potential we assume the M3Y effective interactions [11|. Like in previous 
papers |12|] , we consider only the isoscalar and isovector components for the central 
nucleon-nucleon force : 



vm3y(s = R + r 2 -r 1 ) = v 00 (s) + J 00 5(s) + (v 01 (s) + J i<5(s))ti • r 2 
The isoscalar component of the M3Y force is 



(37) 



voo{r) 



-4r 



7999 



4r 



2134 



-2.5 r 



2.5 r 



MeV 



and the isovector part has the form 
voi(r) = 



-4r -2.5 r 

4885.5 — + 1175.5 



4r 



2.5 r 



MeV 



The strength of the isoscalar zero-range pseudo-exchange potential is taken according 
to the common prescription Joo = -276 MeVfm 3 , whereas the isovector one is Joi = 
217 MeVfm 3 . 

The ground state one-body nuclear densities of the fragments are taken as Fermi 
distributions in the intrinsic frame 



p(r) = po 



1 



1 + exp- [r-Ro 1 + £ PxY X0 {9, 0) 



(38) 



The constant po is fixed by normalizing the proton and neutron density to the Z proton 
and N neutron numbers, respectively. This condition ensures the volume conservation. 
The radius Rq and diffusivity parameters were taken from liquid drop model calculations 
1 13] for the heavy fragments, whereas for the light cluster we consider the prescription 



1/3 

Rq = 1.04 • for the radius and a = 0.35 for the diffusivity. As static deformations, 
[3\, we considered quadrupole, octupole and hexadecupole deformations. 

However, due to the lack of an explicit density dependence in the M3Y effective 
interaction, this potential is characterized by a strong, unphysical attraction of a few 
thousands of MeV for total overlapping. Therefore in the region of nuclear-density over- 



lap we follow the suggestion from [14, 15 1 and introduce a phenomenological repulsive 
potential which originates from the compression effects of the overlapping density 



V Tep {R) = V p J dri J dr 2 pi{ri)p 2 {r 2 )5(s) 



(39) 



where the tildes on the densities signify a distribution of the same shape as (38) but 
possessing an almost sharp surface. The strength of the compression term V p was 
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determined from the nuclear equation of state by requiring for total overlap of two 
nuclei a double normal density of nuclear matter. For a given dinuclear subsystem 
( -Alight + Aheavy) we take the value of the nuclear compression modulus K according 
to the receipt proposed in [16], and compute V p from the equation giving the binding 
energy loss for total overlap |l7]], i.e. R = 0: 

%sy(0) + K cp (0) « J # Ai^ (40) 

The double folded potential (36) is computed by making a general multipole expan- 
sion for two final nuclei with distance R between their center of masses and orientation 
in space given through the Euler angles Qi and f?2 1 15] 



Details on the calculation of the potential multipolar components V^ l ^\ (R) were given 
in§. 

In a previous paper || we considered that when the nuclear molecule is bent the 
reciprocal distances between the heavy fragments and the light cluster are preserved, 
i.e. we allowed the trinuclear molecule to perform only vibrations which result in the 
increase of the angle between the two valence bonds. In this way we excluded possible 
bond stretching vibrations. This time we take into account such a degree of freedom 
and by 5r 13 ( 2 3) we denote the change in the distance between clusters 1(2) and 3. The 
geometry of the system, presented in Fig. 5 for a deformed light cluster, provides the 
values of these quantities as a function of the deflection angles 971,2,3 : 

1 R1R3 , , 2 

Sru = "I(7^my ( '' 1 + ' ^,3, 

*» = -\wrk^-^ (42 > 

On the other hand, from the inspection of Fig. 5, the interaction (|4l| ) between the de- 
formed light cluster 3 and the heavier deformed nucleus 1, reads : 

V(r 13 )= £ J" =^^°(fl 1 +^+gr M )y AlM ( W - e ,0)y Aa> , M ( e +y a ,0) 

XlMM V ( 2A 1 + !)( 2A 2 + 1) 

(43) 

Since we made earlier the approximation ip\ ~ e, the variation of the valence bond 
length ( f42| ) reads 

1 -R1-R3 / Ri + R2 + 2i?3 \ 2 2 fAA \ 
5ri 3 = ~~ — — : — 5 5 : 5 £ (44) 



Expanding the potential (|43|) with respect to the small angle e we have that 

1 



V(r 13 ) = V(R 1 + R 3 ) + ^-Clh 2 (45) 



where 



ri3=#l+#3 

(46) 



12 



Using similar arguments we get the expression for the stiffness of the butterfly mode, 
coming from the interaction between clusters 2 and 3 

23 _ If R \ 2 / R 2 R 3 dV^XM 1 X a xiwooo (r ,\ 

V/ A/ Xi \ z d Ai / r 23 =R2+R3 

(47) 

The last contribution to the stiffness coefficient of the butterfly mode comes from the 
interaction of the heavier fragments 1 and 2. Using again the geometry of Fig. 5 we write 
in multipolar form the interaction between these nuclei 



v{ri - x £ TtgrrW^ ^* 1 + °«*->"°> (48) 



4-7T 

. /7oT- _l 

According to eqs.(|6||8|) the shift in the interfragment distance reads 

Consequently expanding in Taylor series this potential too, we obtain the stiffness coef- 
ficient 

C £ = \ £ fA 1 (A 1 + l)|^-A 2 (A 2 + l) (i?1 - i?2)(i?1 + i?3) 



Ai A2A3 

•0 11 



A3(A3 + l)|^jU A ^ 3 (i?o) 
^1 + ^3 , R R 3 (RiR 3 + R 2 R 3 + 2R 1 R 2 )\ 9V^$JM 



Ai A2A3 



U ° ^ \R 2 + Rz ^R 1+ R 3 ){R 2 + R 3 f J dr lDUj 



In order to obtain the stiffness coefficients of the bond streching vibrations C\ 3 and 
C23, we expand the potentials V{r\ 3 ) and ^(^23) , up to second power of 5r\ 3 and 8r 23 . 
Such an expansion is possible in view of the relative minimum in the potential with 
respect to the inter-cluster distance. As an example we give in Fig. 6 the potential 
between the heavy cluster 1(2) and the light cluster 3. After some algebra we obtain for 
the dinuclear subensemble (13) 

V{r 13 ) = V(Ri + R 3 ) + A 13 5r 13 + ^C 13 5r 2 13 (51) 



where 



/ A 2 (A 2 + 1) R,R 3 dV^ x X(m) \ 

13 " .f, I — 2 — i?TTi^ FAlA2A3(ri3) + — d^ 3 — J _ (52) 

Ai,A 2 ,A 3 V 1,3 / ri 3 =-Ri+R 3 

+ tr + i)^^#^ + ^iLM) (53) 



-R1-R3 9ri 3 <9r 2 
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ri3=-Ri+-R 3 



Since the linear term in Sri 3 has only the effect to shift the origin of the harmonic 
oscillator, the stiffness coefficient of the bond-stretching vibrations is specified by C\ 3 . 
In the same manner we derive the coefficient C23 and next using the relations from 
eq.(|i~8l) we derive the stiffness coefficients of the r and £-modes. 
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4. Applications to three-nuclear molecules 



The model is applied to the systems 96 Sr + 10 Be + 146 Ba, 90 Y + 10 Be + 142 Cs and 
108 Mo + 10 Be + 134 Te. For the computation of the numbers we use me 2 = 938 MeV 
for the nuclear mass, he = 197.33 MeV fm 2 , 1.2 A 1 / 3 for the spherical equivalent radius 
of a heavy cluster and 1.3 A 1 / 3 for the light one. (Because we do not use here a diffuse 
surface as in the last section, for the spherical equivalent radii the usual formula vqAs 
with tq = 1.2 for light and 1.3 for heavy nuclei is used.) In order to obtain 3 j- we set 
it equal to the energy of the 2f state of the individual clusters. 

In case we consider a spherical nucleus, the E(2f) value is taken as the vibrational 
energy. This has to be used with caution because for a rotor the 2f state is a rotational 
one. In a nuclear molecule the rotation of individual clusters is constrained due to the 
link to other clusters. Its rotation is converted into the butterfly motion. This is not the 
case for a vibrational state. In this contribution we will assume a deformed 10 Be where 
the deformation is not taken from Ref. [19] because the assumptions used there are no 
longer valid for light deformed nuclei. We rather use the SU(3) model of the nucleus 
and deduced from there a deformation of 0.175 (see Ref. 0]). It will be seen that the 
results do not sensitively depend on that. Interpreting the 2f state at 3.368 MeV in 
10 Be as rotational, it is absorbed into the butterfly motion. Instead the first vibrational 
state is a 7-mode at 5.958 MeV, indicating a very stiff system. The fact that the 3.368 
MeV transition (minus the 6 keV shift) is seen in experiment [l], Q] speaks in favour of 
a vibrational 10 Be nucleus. We carried out computations in frame of the Hartree-Fock 
method, with pairing correlations taken into account and using Skyrme III forces. The 
result was that the deformation energy curve of 10 Be has a spherical minimium and it 
is symmetric for small deformations. Nevertheless, we would like to see the effects of a 
possibly deformed light cluster. All this has to be taken into account when it comes to 
the interpretation of the theoretical results. 

After having described the three systems mentioned above, shortly the nuclear struc- 
ture of the participants of the system 90 Y + 10 Be + 142 Cs will be discussed. 90 Y and 
142 Cs are odd-even nuclei with an odd number of protons. Their treatment requires the 
inclusion of the spins of the extra protons. Though, the extension can be done it would 
break the scope of this contribution. In a subsequent paper it will be treated. The com- 
ments we will make, however, are important for the interpretation of the experimental 
results. 



4.1 96 Sr + 10 Be + 146 Be 

The nuclear structure of this system was already discussed in Refs. M, Hj. The new 
contribution here is the inclusion of the (3 and 7 and the relative vibrational modes. 
Although, the relative vibrational modes where discussed in Ref. Q numerical values 
where not given due to the missing information on the stiffness parameters. These will 
be provided now. 

The Sr and the Ba nuclei are prolately deformed and the corresponding deformation 
values are given in Table 1. As outlined in the last section and using the data from Ref. 
@ the parameters of the Hamiltonian are deduced and listed in Table 1. The spectrum 
can then be determined from the Eq. (|34]). The radii along the prolate symmetry axis 
are also listed in Table 1. For the Be nucleus we showed in Ref. pj that according 
to the SU{3) model it can be taken as triaxial and the deformation is 0.175. This 
consideration does exclude any mixing with other SU(3) representations due to the 
SU(3) mixing terms like pairing and spin-orbit interaction. There are many indications 
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that the 10 Be nucleus can be assumed to be spherical Q. However, in order to see the 
influence of a deformed light cluster we will assume a deformed 10 Be nucleus. 

Using the parameters of this system listed in Table 1, the spectrum of the molecular 
states is given by (units are in MeV) 



E = 0.000862[/(J + 1) - K 2 } + 2.406(| K-Kx-K 2 -K 3 \ +2n e ) 
+1.229710! + 1.05371/32 + 6.179n A 

+1.507(- | K x | +2n m ) + 1.566(- | K 2 \ +2n m ) + 5.958(^ | K 3 \ +2n m ) 
+3.61n r + 17.59n ? , (54) 

where we have skipped the zero point energy contribution, i.e. E gives the difference 
in energy to the ground state. All deformation vibrational states are lying above 1 
MeV. The same holds for the butterfly frequency. In conclusion, below 1 MeV only the 
rotational states belonging to the ground state band appear. The relative £ motion is 
at such a large energy that it does not play any practical role. Interesting to note is 
that in the calculation of the dominant contribution comes from the last term of 
011 as given in Eq. (|lq). The other terms contribute at most three percent. Even the 
rotational contribution from the light Be cluster, given by the term before the last one 
of 0ii, can be neglected due to a small £o- This is similar to the two cluster case where 
the corresponding term is dominating all others. 

In Fig. 7 the expected structure of the spectrum is plotted. Only band heads con- 



tained in Eq. (54) are shown and not those belonging to other degrees of freedom, like 
the rotational octupole band head state 1~ in 146 Ba and the extra + state at 1.465 
MeV in 96 Sr. As can be seen in Fig. 7, there are no excited band head states below 1 
MeV suggesting a stable behavior against the butterfly motion. Nearly all states below 
1 MeV are rotational one belonging to the ground state band. Note that in Ref. [|j The 
hu> £ was estimated assuming a spherical 10 Be nucleus. 

The energy values of the first 2f and 4f states are 5.2 keV and 17 keV respectively. 
The deformation of 10 Be results in an increase of the separation of the heavy clusters 
which raises the moment of inertia. The expected lowering in the energies of the rota- 
tional states is small compared to the results of Ref. Q, which are 6 keV and 20 keV 
for the 2^ and 4^" states, indicating a small influence of the supposed deformation of 
the 10 Be nucleus. 



4.2 112 Ru + 10 Be + 130 Sn 

This system is the most symmetric one we could get for which experimental information 
about the structure of the individual clusters are available and not just the ground state 
only. This is important for deducing the deformation of the nuclei. This system has not 
been seen yet but should exist. 

The heavy fragments are again even-even nuclei. Using the tables of Ref. [|l9[ the 
corrected deformation of 112 Ru is given by 0.237 corresponding to a large j3. For 130 Sn 
no information in these tables are available. However, The Sn isotopes are known to be 
an excellent example for the seniority scheme The proton shell is closed and the 
neutron shell is open. Because the seniority scheme is realized a zero deformation can 



be assumed. This is also confirmed by the {E{^[) / E{2^)) ratio [22] which is 1.63 for 
130 Sn. For 112 Ru the ratio is 2.72 indicating a rotational structure. 
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The parameters of the nuclei and the system are listed in Table 1, including the radii 
of the clusters along the line of contacts. The excitation energy is given by (units are 
in MeV) 

E = 0.00094[J(J + 1) - K 2 } + 2.215(| K-K 1 -K 2 -K 3 \ +2n £ ) 
+0.524(i | K\ | +2n Vl ) + 5.958(^ | K 3 \ +2n m ) 
+3.911n r + 19.035n 5 

+1.220iV Sn , (55) 

where the last term describes the five-dimensional harmonic oscillator for the Sn nucleus. 
The spectrum is presented in Fig. 8, with the same characteristics as in 96 Sr + 10 Be + 
146 Ba. Not plotted are the band heads belonging to different degrees of freedom than 
those described by the model (see discussion in the former section). Also suppressed are 
the states of the five dimensional harmonic oscillator of 130 Sn with Tuop = 1.220 MeV. 
There is a 7 vibrational state at approximately half a MeV though. One should observe 
above this band head a rotational structure similar to the ground state band with the 
difference that also a 3 + state exists at about 5.5 keV and a 4 + state at 12.7 keV above 
the 2 + band head state. The £ relative vibrational term can also be neglected. In this 
system the rotational part is dominated to almost 100% by the last term of On of Eq. 
(|l6|). The influence of the assumed deformation of 10 Be is again small. 



4.3 108 Mo + 10 Be + 134 Te 

The heavy fragments are even-even nuclei. Using the tables of Ref. the deformation 
of 108 Mo is given by 0.354 which we corrected to 0.333 using the additional deformation 
dependent terms in the B(E2;0f — > 2^) as given in Ref. [pXj ] . The experimental 



information of the spectrum is taken from Ref. [20|. For 134 Te no information in these 
tables are available, however, the tendency observed coming from the lighter isotopes 
indicates a small f3. When we look at the ratio (E(4f)/E(2f)) we obtain 2.92 and 1.23 
for 108 Mo and 134 Te respectively. This experimental observation supports a deformed 
nucleus for 108 Mo and possibly a spherical deformation for 134 Te. 

The 108 Mo nucleus is particularly difficult to treat. Within the geometrical model 
1C] the Potential-Energy-Surface (PES) of the neigbouring nucleus 108 Ru has a spherical 



absolute minimum and a triaxial local minimum at large deformation 23]. The energy 
of the ground state, however, lies above the saddle point and the ground state is a strong 
mixture between both deformations. A large (3 in average is also indicated by the results 
of (3f3 decay for 108 Mo where with great success a model for strongly deformed nuclei 



was applied |24| . 

Under the assumption that the Mo isotope is deformed while the Te nucleus is a 
vibrator and using the parameters listed in Table 1, the excitation energy is given by 



E = 0.00092[/(I + 1)-K 2 } + 2.127(| K-K x - K 2 - K 3 \ +2n e ) 

+0.586(i I K x I +2n m ) + 5.958(^ | K 3 | +2n m ) 
+3.54n r + 17.30n ? 

+1.280V Te , (56) 

where the last term describes the five-dimensional harmonic oscillator for the Te nucleus. 
The spectrum is presented in Fig. 9, again with the same characteristics as in 96 Sr 
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+ Be + Ba. Not plotted are the band heads belonging to degrees of freedom 
different from those of the model. Also excluded are the states of the five dimensional 
harmonic oscillator of 134 Te with FnuTe = 1-280 MeV. There is a 7 vibrational state at 
approximately half a MeV though. One should observe above this band head a rotational 
structure similar to the ground state band with the difference that also a 3 + state exists 
at about 5.5keV and a 4 + state at 12.8keV above the 2 + band head state. The £ relative 
vibrational term can also be neglected. In this system the rotational part is dominated 
to almost 100% by the last term of 6n of Eq. (|l6|). Again the interpretation with 
respect to the influence of the 10 Be deformations similar as in the two former cases. 



4.4 Structure of the participants in 

90 Y + 10 Be + 142 Cs 

This is the third system possibly identified in the experiment of Ref. ||. The heavy 
clusters are odd-odd nuclei. The influence of the extra odd protons should be included 
in our consideration. The way to do it for individual nuclei is given in Ref. Jl^]. It 
would, however, exceed the scope of this contribution to add a corresponding discussion. 
We refer to a later publication on this subject. Therefore, we only discuss the nuclear 
structure of the participants. We consider it important because the interpretation of 
the experimental results rely heavily on the structure of the heavy clusters. They can 
be described by an even-even core and an odd proton around it |E0"|| . The deformation 
can be deduced via the one of neighboring nuclei. Take -jgS^, toZrso as neighboring 
nuclei and also IgSrso as the core. The deformation of these nuclei are listed in Ref. ||, 
except for ggS^. The deformation are 0.09 and 0.12 for the last two nuclei respectively. 
However, one must take some care in using these tables. The formula used to deduce the 
deformation from the experimental B(E2,0f — > 2^) transition is theoretically biased, 
i.e. a strongly axial symmetric deformation is assumed. In this case the deformation 
value is proportional to the square root of the transition. For spherical and triaxial 
nuclei the formula is incorrect and in general the deformations deduced are too high due 
to the fact that higher orders ]l~0[ | in deformation are also neglected. Nevertheless, the 
tables in Ref. jl{| give a good idea about the trends. Another possibility is to use the 
tables of M. Sakai [22] where the ratio (E(4f)/E(2^)) is investigated. For rotational 



nuclei the ratio is 3.33 and for vibrational nuclei it is 2.0. Most nuclei lie in between 
these two values. For the three nuclei mentioned above this ratio is 1.99, 1.41 and 3.15 
respectively. Except for the last value, which is near to a rotator, the data indicate a 
spherical deformation, in agreement with the data listed in Ref. [19]. 

For 5| 2 Cs87 as neighboring nuclei we took g| 2 Xegg, 5g 2 Bag6 and for the core 54°Xe86- 



The deformation values listed in Ref. [19| are 0.157 and 0.114 for the last two systems. 
No information is listed for the first nucleus. Using the tables of Sakai |[22|| , except for the 
first nucleus were we used the ISOTOPE EXPLORER @, the ratio of (E(4f)/E(2^)) 
are respectively 2.41, 2.32 and 2.22. They hint to a spherical nucleus while the (3 values 
are in between. In Ref. || a deformed nucleus was assumed, which are not confirmed 
by data. 



5. Conclusions 

In this contribution we proposed a model for three cluster molecules, whose existence are 
suggested by experiments in Refs. ]l], ||. First results of their theoretical description are 
published in Refs. ||, ^]. The modes of this molecule are the butterfly mode, rotations 
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of the whole system and (3 and 7 vibrations of each cluster. The belly-dancer mode is a 
particular rotation around the molecular z-axis. In general, all clusters can be deformed. 

In order to render the model soluble, assumptions had to be introduced. First, the 
three clusters are supposed to form linear chain molecule with the lightest cluster in the 
middle. Second, the inclination of the symmetry axis of the clusters and the distance 
of the lightest cluster to the molecular z axis were assumed to be small. For later 
applications we restrict the consideration to the case where the first cluster touches the 
third one and this the second nucleus. Relative vibrations were also added but taken 
to be in the direction of the z axis. Even with these assumptions, for large asymmetric 
systems there is a coupling between the butterfly motion (denoted by the angle e) and 
the rotation around the molecular y axis, which is perpendicular to the (x, z) plane in 
which the center of masses of the three clusters are located. For "nearly" symmetric 
clusters and when the ratio mi ™ 3 m2 is very small, this coupling can also be neglected. 
When there is a stretching motion along the molecular x axis, coupling terms appear 
which have to be diagonalized in the basis constructed in this contribution. 

One basic feature of the model is that the butterfly motion is alway possible and 
the contributions come from the motion of the center of masses of the three clusters, 
contrary to the two cluster case where the deformation plays the essential role for Tilo £ . 
The contribution to is more than half from the center of mass terms in Eq. (16) and 
less than half coming from the deformation dependent terms. In the last terms the 10 Be 
contribution is negligible. 

The model was applied to the systems 96 Sr + 10 Be + 146 Ba, 122 Sn + 10 Be + 120 Ru, 
and 108 Mo + 10 Be + 134 Te. All these splittings have the characteristic that there 
are few vibrational states below 1 MeV. The last two systems exhibit a low lying 7 
vibrational band head at about 0.5 MeV. The rotational states of the ground state band 
are strongly squeezed. The first excited 2 + state is at approximately 5-6 keV and the 
4|" state at approximately 17-19 keV. The butterfly motion lies at energies larger than 
3 MeV indicating a particular strong stiffness against this mode, contrary to the two 
cluster molecule. The structure, obtained for these systems, have to be considered as 
a rough approximation, due to the assumptions made and the neglection of couplings 
between most of the degrees of freedoms. 

The structure obtained for individual nuclei of the last system and for 90 Y + 10 Be + 
142 Cs, where only a discussion on the nuclear structure of the individual clusters is given, 
are in contradiction to Ref. Q. There, for example the 142 Cs nucleus is assumed to be 
deformed while experimental data prefer a spherical nucleus in its ground state. Also 
the nucleus 108 Mo is taken as spherical in Ref. while the data suggest a deformed 
nucleus. Because for the three cluster molecule the main contributions come from their 
center of masses, changes in the deformation of the nuclei will not affect the main results. 
However, in Ref. ||] the deformation plays an important role because it explains the 
possible systematics observed of the shift in the energy of the 2f state in 10 Be. Based 
on the structural investigation we made, the explanation of the shift is probably not as 
easy. We suggest explicit microscopic investigations of the influence of a heavy cluster 
to the smaller one, maybe implying a change in structure of also the heavy fragments 
from deformed to spherical and vice versa. Also a revision of the experimental analysis 
is recommended. 
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APPENDIX A: 

The scalar product (^ lab • ^ lab ) in terms of the ve- 
locities uj[ and i: 

As spherical components we use the expression 



1 

£o = iz (57) 



where the definition of £±i differs from the usual one || for convenience. In the molecular 
frame the vector £ is defined to lie in the molecular plane given by the molecular z axis, 
parallel to the vector r. Therefore, the component vanishes and the relation of the 
cartesian to the spherical components is such that the £±i components are in their 
absolute value identical and are given by ±-^=£ x . 

Neglecting terns of the order 0(e 3 ) in Eq.(9) we obtain 

& = Co+ 2 R 2 + R 3 A 1 + A 2 £ (58) 

where £0 is the value of the ^-coordinate in the linear chain configuration (see eq.(10)). 
The time derivative of is determined using the procedure outlined in Ref . ju]] . After 
some algebra, which implies also the calculation of Wigner D-functions time derivatives 
we arrive at 



lab -i lab ) = J2(- 1 ) m ^- h m « (Ri + R 3 ) 2 e 2 + e (u;' 1 2 +u;'i) + (R 1 +R 3 ) 2 s 2 u' 3 2 

m 

+ 2^{R l + R z ){iu' 2 -£J 1 J 2 ) . (59) 

APPENDIX B: 

Coordinate symmetries of the tri-nuclear 
molecule with prolate nuclei 

We follow the procedure outlined in the Refs. |7|, 10]. The transformations of the 
molecular frame consists of the operators R\ m and R\ ml the m refers to the molecular 
frame, and their action is given by 

Ri,m(x,y,z) = (x,-y,-z) 
R2,m(x,y,z) = (-x,-y,z) . 

(60) 
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These operators act on the Euler angles (k = 1,2,3), the other coordinates \i-> <\>ii 
<3?i {% = 1,2,3) and f and £. The result is given in Table 2. The angles Xii 4>i an d &i 
correspond to the Euler angles describing the rotation of the nucleus from the molecular 
frame to the principal axis of cluster no. i. 

Because the three clusters are supposed to lie in a plane the angles Xi are P ut to 
zero. Also a small butterfly angle is assumed and all angles ipi are proportional to it. A 
coordinate symmetry transformation consists of those actions where the angels Xi an d 
ifi are changed at most by a sign. Inspecting Table 2, the only combinations allowed 
are Ri >m Ri >Pi and i?2p» whose action is given in Table 3. 

These operators have to be applied to the solution 

M?i,02,03,e,r,fr,ft,&,f7i,^ 

X9n ( (O9n 01 (Pl)g n02 (h)9n H {h)XK x ,n m (Vl)XK 2 ,n V2 (m)XK 2 ,n V2 (m) ( 61 ) 

of the Schrodinger equation. The action of R.2, Pi (i = 1, 2, 3) on this state changes 3>j to 
$j ± | and rji to —rji, i.e. it acts only on e tKl<i ' i and XKi^XVi)- The result is a phase 
(— l) Ki implying only even values of Ki. 

The action of the operator Ri jjn Ri, Pi is more involved. The result is 

N $2, $3) + (-ly-XDif.KWfi-Qu -$ 2 , -<& 3 )} 

XXK,ne( £ )9n r (f)g nii {€)9n [h (Pl)9n 02 {h)9n 03 (Jk)XK unrn (Vl)XK 2 ,n V2 {m)XK 2 ,n V2 

(62) 

where J\f is a normalization factor and 

/(<&! $2 $3) = e^^i^i+^^a+^^a) _|_ e -i{Ki^ 1 +K 2 'i>2-K 3 ^ s ) 

^(K^l- K 2 § 2 - Kz$z) _|_ e -j(-ft"l*l-A"2"I > 2+-ft"3 < I > 3) (g3) 

The quantum numbers acquire the possible values 



Ki = 0,2,4,... 

K = 0,1,2,3,4,... 

L = K,K + l,K + 2,... 

n r , n ? , n £ , n m = 0,1,2,3,4,... . (64) 
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FIGURE CAPTIONS: 



Figure 1: Illustration of the main variables of a three cluster molecule. The light cluster 
is plotted as spherical. In general it will be deformed as the heavy nuclei. 

Figure 2: In the upper half the mode discussed in this paper is plotted. The lower half 
presents a more complex case which corresponds to the anti butterfly mode in the two 
nuclear molecule for the limit of vanishing mass of the light cluster. 

Figure 3: The potential of the light cluster 3 in the field of the two heavier fragments 
along the molecular z axis, for three fixed tip distances : d= 2 (solid lines), 3 (dashed 
line), 4 (dotted line) FM. The trinuclear molecule comprises 96 Sr + 10 Be + 146 Ba. 

Figure 4: Main variables for the discussion of the motion of the clusters with respect to 
each other. Only the centers of mass of the nuclei are plotted. 

Figure 5: In case the third cluster is deformed, the alignment of a heavy nucleus with 
the light one is not perfect. The relations of the angles are illustrated in this figure. 

Figure 6: The potential between the heavy fragment 146 Ba and 10 Be (solid line) and 
between 96 Sr and 10 Be (dashed line). 

Figure 7: Spectrum of the system 96 Sr + 10 Be + 146 Ba. For detailed explanations, see 
the text. Only band heads are shown. On top of each band head there is a rotational 
band with the characteristics explained in the text. The butterfly mode is the 1 + state 
to the right and the first relative vibration is given by the 1" state. The 2 + states 
belong to the 7 vibrational states with either K\ or Ki equal to 2. The + band heads 
consist of the ground state, (3 and 7 vibrational (with n Vi = 1) band heads of the heavy 
clusters. 

Figure 8: Spectrum of the system 112 Ru + 10 Be + 130 Sn. On top of each band head 
there is a rotational band with the characteristics explained in the text. The butterfly 
mode is the 1 + state to the right and the first relative vibration is given by the 1~ state. 
The 2 + state belongs to the 7 vibrational state of the first cluster with K\ = 2. The 
+ band heads consist of the ground state and the 7 vibrational (with n m = 1,...,4) 
band heads of 112 Ru. No (3 vibrational states could be identified in the three clusters at 
energies below 5 MeV. 

Figure 9: Spectrum of the system 108 Mo + 10 Be + 134 Te. On top of each band head 
there is a rotational band with the characteristics explained in the text. The butterfly 
mode is the 1 + state to the right and the first relative vibration is given by the 1~ state. 
The 2 + state belongs to the 7 vibrational state of the first cluster with K\ = 2. The 
+ band heads consist of the ground state and the 7 vibrational (with n m = 1,...,3) 
band heads of 108 Mo. No f3 vibrational states could be identified in the three clusters 
at energies below 5 MeV. 
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system 


™Sr + 1U Be + u « Ba 


112 Ru + 1U Be Sn 


IU8 Mo + I0 £> e + 134 Te 


An 


0.338 


0.237 


0.333 


A)2 


0.199 


0. 


0. 


A)3 

/ uo 


0.175 


0.175 


0.175 


a 2 

2011 


8.62 10 -4 


9.10 10~ 4 


9.20 10~ 4 


/ LUJf 


4 970 


3 Q1 1 


3 850 




1 Q 41 7 




1 8 ^1 1 




9 497 




9 

Zi . UO J- 


't"^/3i 


1 99Q 






TiUln 


1.053 


1.220 


1.280 




6.179 


6.179 


6.179 


hu! rn 


1.507 


0.524 


0.586 




1.566 






hoo m 


5.958 


5.958 


5.958 


Ri 


6.66 


6.88 


6.92 


R2 


7.11 


6.08 


6.14 


R 3 


3.11 


3.11 


3.11 



Table 1: Parameters of the three systems discussed in this paper. In case the deformation (3k, 
for one particular nucleus k, is zero the Tiujp k has to be interpreted as the energy hu k of the 
five dimensional harmonic oscillator. For the case of an oscillator the corresponding hu Vk is 
put to zero because it is not relevant. The deformation parameters have no units. The one 
of the hujk are in MeV. The units of the radii are is in fm, where we used for the spherical 
equivalent radius the formula r^A 1 ^ with r = 1.2 for a heavy and r = 1.3 for the light 
cluster. An "-" indicates that either no information is available, very insecure or /3ok is zero. 
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z 


z 


Z 
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Table 2: The action of the basic coordinate symmetry operators on the collective variables 
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Table 3: The action of the allowed combinations of symmetry operators which satisfy the 
condition that after their application \i is still zero and ifi is maintained near zero. 
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